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ABSTRACT 

GRB 070125 is among the most energetic bursts detected and the most extensively 
observed so far. Nevertheless, unresolved issues are still open in the literature on the 
physics of the afterglow and on the GRB environment. In particular. GRB 070125 
was claimed to have exploded in a galactic halo environment, based on the unique- 
ness of the optical spectrum and the non-detection of an underlying host galaxy. In 
this work we collect all publicly available data and address these issues by modelling 
the NIR-to-X-ray spectral energy distribution (SED) and studying the high signal- 
to-noise VLT/FORS afterglow spectrum in comparison with a larger sample of GRB 
absorbers. The SED reveals a synchrotron cooling break in the UV, low equivalent 
hydrogen column density and little reddening caused by a LMC- or SMC-type extinc- 
tion curve. From the weak Mgn absorption at z = 1.5477 in the spectrum, we derived 
log7V(Mgii) = 12.96lg'}g and upper limits on the ionic column density of several 
metals. These suggest that the GRB absorber is most likely a Lyman limit system 
with a 0.03 Z < Z < 1.3 Z & metallicity. The comparison with other GRB absorbers 
places GRB 070125 at the low end of the absorption line equivalent width distribution, 
confirming that weak spectral features and spectral-line poor absorbers are not so un- 
common in afterglow spectra. Moreover, we show that the effect of photo-ionization 
on the gas surrounding the GRB, combined with a low iV(Hi) along a short segment 
of the line of sight within the host galaxy, can explain the lack of spectral features in 
GRB 070125. Finally, the non-detection of an underlying galaxy is consistent with a 
faint GRB host galaxy, well within the GRB host brightness distribution. Thus, the 
possibility that GRB 070125 is simply located in the outskirts of a gas-rich, massive 
star-forming region inside its small and faint host galaxy seems more likely than a 
gas-poor, halo environment origin. 

Key words: gamma-ray burst: individual: GRB 070125 



1 INTRODUCTION 

The progenitors of long gamma ray bursts (GRBs) are 

* Based on target-of-opportunity observations collected in ser- belie ved to be rapid ly rotating massive stars ( for a review, 

vice mode under program ID 078.D-0416, P.I. Vreeswijk, with see IWooslevI ll993J ; iMacFadven fc Wooslevl 119991 ). The 
the FOcal Reducer /lo w dispersion Spectrograph 2 (FORS2 
lAppenzeller et al,lll99Sr) installed at the Cassegrain focus of the 

Very Large Telescope (VLT), Unit 1, Antu, operated by the Eu- For further information or questions on the content of the paper, 

ropean Southern Observatory (ESO) on Cerro Paranal in Chile. please e-mail annalisa@raunvis.hi.is. 
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observations of some SNe asso ciated with GRB s have 
corroborated this hy pothesis IjWooslev fc Blooml 120061 ; 
iHiorth fc Blooml l201ll . and references therein), implying 
that long GRBs are loc ated in star-f o rming reg i ons. W ithin 
the fireball model (see IPiranl |2004| ; iMeszarosi 120061 , for a 
review) the GRB afterglow is produced by the interaction 
of an expanding blast wave with the surrounding medium. 
The external shocks are responsible for the synchrotron 
radiation, from X-rays to optical and down to radio fre- 
quencies as the blast wave decelerates. The spectral energy 
distribution (S EP) and its evolution is described by the 
fireba ll model (pari. Piran fc Naravanlll998l ; iGranot fc Saril 
2002). Modelling the SED can simultaneously provide 
information on the intrinsic spectral shape of the GRB, i.e. 
the fireball physical parameters, and on the GRB local and 
host galaxy environment. In the latter case, this includes 
the circumburst density or the line of sight dust extinction, 
depending on whether the broadban d blast wave, down to 
radio fr equencies, is modelled (e.g., Panaitescu fc Kumarl 
2002; Biornsson. Gudmundsson fc Jo hanncssonl 12004 



dc Ugartc Postigo ct al. 2005: 



van dor H orst et al.l 2005; 
Johannesson. Bi ornsson & Gudmundsson 2006; Rol et al. 
2007|) or the optical-to - X-ray SEP llStarling et al.l 12001 
20081 : ISchadv et alJlJOlcl ; IZafar et al.ll201lh . The interstellar 



medium of the GRB host galaxy can also be investigated 
through absorption line spectroscopy of the afterglow, 
providing the chemical composition of different ions 
the absorbing gas and 



Vreeswiik et al.l 2004 ; iFvnbo 



ll 



Prochaska et al.ll2007l ; iLedoux et alj|2009r i 



possibly its metallicity (e.g. , 
jo et all 120061 : ISavagliol 120061; 



GRB 070125 is among the most energetic bursts de- 
tected so far, both in the prompt high energy release (with 
an is otropic emitted energy i?i ao ~ 10 54 erg: iBellm et al] 
12008! ) and its afterglow (with a collimation- corrected blast 
wave kinetic energy of more than 10 52 erg: IChandra et al] 
|2008i ). The optical afterglow is amongst the brightest at 
around one day after the burst, when redshift-corrected 
(R c ~ 17 shifted at z = 1: lUpdike et al"1l200Sl ; iKann et all 
|2010j). The brightness of the event allowed a global mon- 
itoring campaign of the afterglow resulting in one of the 
largest datasets collected so far, covering radio, mm, NIR, 
optical, UV and X-rays ( Cenko ct al. 2008; Chandr a et al] 



l200Sl : iPai et alj|2008l : lUpdike et alj|2008l : IKann et al.ll2010l ) . 
Pespite the observational effort, there is still much debate 
over both the jet properties and the type of environment in 
which GRB 070125 was situated. 

Based on the absence of strong spectral features and of a 
bright underlying host, GRB 070125 was reported to be the 
first long GRB to have exploded in a galactic halo environ- 
ment. This would pose a challenge to modelling of progenitor 
objects. In particular, the GRB site was suggested to be in 
a compact stellar cluster that resulting from the interaction 
of two blue fai nt galaxies at a pr ojected distance of ~ 27 
and « 46 kpqj (|Cenko et alJbO OS'). However, a very limited 
sample of GRB afterglows was used as reference to argue 
the peculiarity of the spectrum (few spectral features, weak 
Mg 11 absorption and low inferred iV( Hl)). By contras t , from 
a partial broadband SEP analysis, IChandra et al] ((2008) 



1 Derived from the 3.2" and 5.5" angular distance between the 
afterglow and the two putative host galaxies. 



found a very high circumburst density (n ~ 50 cm -3 , from 
the kinetic energy and n « 40 cm -3 , from the broadband 
modelling of the SEP). However, in their data it was not 
possible to distinguish between a constant or a wind-like 
density profile nor to firmly constrain the putative (chro- 
matic) jet break an d the inve r se Co mpton influence on it. 
On the other hand, IPai et al] (j2008i ) claimed a "textbook" 
achromatic jet break at 5.8 days after the burst. Another 
open issue is the presence or absence of a synchrotron cool- 
ing break in the SEP between the optical and X-rays and 
the agreement or not of the data w ith the fireball closure, 
as discussed in (|Updike et all 2008'). 

In this work, we collect the multitude of observations 
spread through the literature and analyse them together in 
order to obtain a coherent picture of the GRB environment. 
In particular, we investigate the location of a synchrotron 
cooling break and the environment dust and metal content 
by studying the NIR-to-X-ray SEP at different epochs. Fur- 
thermore, we analyse the very high signal-to-noise (S/N) 
VLT/FORS spectrum of the afterglow and discuss the possi- 
ble implications for the absorber environment. We then com- 
pare the spectral features to th e large sampl e of G RB after- 
glows with redshift collected in IFvnbo et al] (|2009r ) and two 
sub-samples of GRBs potentially similar to GRB 070125. 
The observations and data analysis are presented in Sec. 
[2] and Sec. [3] respectively, the SEP modelling in Sec. [4] and 
the discussion in Sec. [5] Our conclusions are summarized in 
Seed 

Throughout the paper we use the convention F v (t) oc 
t~ a v~^ for the flux density of the afterglow, where a is 
the temporal slope and /3 is the sp ectral slope. We refer 
to the Solar abundances measured bv lAsplund et al] (|2009Q 
and adopt cm" 2 , as the linear unit of column densities, N. 
Hereafter, we assume a standard ACPM cosmology with 
H = 70.4 km s" 1 Mpc" 1 , fi M = 0.27 and fi A = 0.73 
|jarosik et al]|201ll ). 



2 OBSERVATIONS AND DATA REDUCTION 

GRB 070125 was localized by the Interplanetary Network 
(Mars Odyssey, S uzaku, INTEGRAL and RHESSI) at 
to =07:20:45 UT (|Hurlev et al.l I2007J ). and subsequently 
observed by Swift. The high energy properties of th e 
GRB prompt emission are discussed in lBellm et al] J2 008). 
The optical afterglow was first localized by ICenko fc Foxl 
(|2007h at the position R.A. = 07 h 51 m 17.75 s and Pec. = 
+31°09'04.2" (J2000) and the host galaxy redshift, z = 
1.547, was determined from a Mgn doublet absorption in a 
Gemini-North/GMOS afterglow spectrum l|Fox et alj|2007f ). 
Thanks to the brightness of the afterglow, a multitude of 
observations in th e optical and NIR could be secured with 
20 telescopes (e.g . , IChandra et"aLll200Sl : lUpdike et ai1l2008l : 
IKann et alll2010h . 



2.1 Swift XRT 

The Swift X-ray telescope fXRT. lBurrows et al.ll2005T ) began 
observing the GRB 070125 field 12.97 hrs after the trigger 
and continued to monitor the decaying afterglow in photon 
counting (PC) mode for 15 orbits or 173 ks of exposure time. 
The last observations ended on 2007 February 12 at 23:50 
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Figure 1. The radio to X-ray afterglow light curves, scaled with the factors reported in the legend for a clearer display. The dashed 
vertical lines mark the SED times. 



UT, 18.69 days after the burst trigger. We obtained the XRT 
event files from the UK Swift Science Data Centre archivqj 



http : //www. swift . ac.uk/swift_portal 



and processed them using the standard Swift XRT pipeline 
software version 0.12.3 and calibration database CALDB-11. 

We extracted X-ray spectra and response files for each 
orbit individually as well as a late-time combined spectrum 
comprising orbit 3 onwards. Spectra were grouped such that 
a minimum of 20 counts fell in each bin and spectral analy- 
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sis was perform ed with the XSPEC spectral fitting package 
l|Arnaudlll996l ), version 12.6.0. 

We extracted a 0.3-10 keV light curve and transformed 
the count rate to flux using 3 conversion factors derived 
from spectra of orbit 1, orbit 2 and orbit 3 onwards of 1 
ct s _1 = 4.44,4.09,3.67 x 10" 11 erg cm~ 2 s" 1 , respectively. 
These three intervals were selected in order to apply the 
proper conversion factor to each regime, from the brighter 
to the late fainter epochs. We also collected the dynami- 
cally binned XRT flux light curve from the Swift reposi- 
tory (jEvans et al.l 120071 ). which was converted from count 
rate to flux with the average conversion factor, 1 ct s _1 
= 3.8 x 10 _11 ergs cm -2 s _1 . We then converted the binned 
and unbinned light curves into flux density F v at 1.73 keV, 
the logarithmic average of the XRT band (0.3-10 keV), as- 
suming a spectral slope /3 = 1.1, as found in Sect. |3~T1 The 
binned XRT light curve is shown in Fig. [T] 



2.2 Swift UVOT 

The afterglow of GRB 070125 was d etected with the Sw ift 
ultraviolet/optical telescope ("LTVOT. lRoming et ajj2005J ) in 
all six lenticular filters spanning v to uvw2 (covering a range 
between 546-193 nm). These observations were coincident 
with those of the Swift XRT during the period 12.97 hrs to 
18.69 days after the trigger. The optical/UV transient was 
detected out to 3.5 days after the trigger. We inspected the 
UVOT sky images and selected those of sufficient quality 
to perform photometry to an accuracy of 0.4 magnitudes or 
better. We performed photometry on the images using the 
UVOTSOURCE task included in HEAS OFT v6.7 and calibra- 
tions provided bv lPoole et al.l 1)20081 ). We used a 5" radius 
source extraction region, and a nearby 20" radius source- 
free background region. The UVOT photometry is reported 
in the appendix 



2.3 VLT/FORS2 

Starting at 04:25 UT on 2007 January 26 (21.1 hours post- 
burst), a series of spectra were obtained with VLT/FORS2 
in long-slit spectroscopy mode with a l'.'O wide slit. The se- 
quence of grisms used was 300V, 600z+OG590, 1400V and 
1200R+GG435. This allowed us to both cover a larger wave- 
length window with the lower resolution grism (300V) and 
obtain mid-resolution spectroscopy for different regions of 
the spectrum. Cosmic ray removal was performed on each 
individual spectrum u sing the Lapla c ian C osmic Ray Iden- 
tification algorithm of Ivan Dokkuml (|2001l ). The seeing re- 
mained relatively stable during the observations, between 
0'.'6 and 079, yielding the spectral resolutions and excellent 
S/N, reported in Tabled] 



2.4 Literature data collection 

We gathered all the data on GRB 70125 available in th e lit- 

erature from 

iKann et al.l I 

exclude the observations from GRB Coordin a tes N etwork 

(GCNs) circulars reported in IChandra et al.l l|2008l ). with 



the exception of three NIR datapointqj (K s , J, H), given 
the limited NIR coverage. In addition, we include the Swift 
UVOT and XRT observations that we analysed. In all, 365 
NIR to UV photometric data points were collected, for a 
total of 19 filters, and 77 radio and mm data points. We 
converted the magnitudes into flux density using the zero- 
points available for each specific filter. All the magnitudes 
were corrected for Galac tic extinction (E(B — V) = 0.052 
mag), as listed bv lSchlegel. Finkbeiner fc Davis! (1 19981 ). and 
for transmission through the Lya for est at redshif t 1.5477 
in each optical and UV band (e.g., iMadaul 1 19951 ) adopt- 
ing the s pectral slope, /3=0 .58, as derived from the opti- 
cal SED (jUpdike et al.ll2008l ). The corresponding transmis- 
sions are tru = 0.996, tr uvw i = 0.848, tr uvm 2 = 0.633 and 
tr UVW 2 — 0.539. The light curves for each filter are shown in 
Fig-IU The dataset that we used is rather inhomogeneous, as 
it is collected from different sources in the literature. Despite 
the uncertainty on the single data points, the interpolation 
of the flux densities at a given epoch are weighted by these 
errors, giving more importance to the more precise measure- 
ments. Thus, the uncertainties due to the inhomogeneous 
sample are limited, given the large dataset, and cannot sig- 
nificantly influence our results. 



3 SPECTRAL ANALYSIS 

3.1 X-ray spectrum 

We first independently modelled the spectrum of each or- 
bit with a simple power law, absorbed by gas at redshift 
Z = 1.5477, as derived from the spectral analysis (see Sec. 
13. 2|) . finding no evidence for spectral variability in the XRT 
data alone within this mod el. The Galactic abso rption was 
fixed at 4.3 x 10 20 cm -2 JKalberla et al.ll2005l ). We then 
extracted the time-averaged spectrum over the whole XRT 
observation. Modelling the overall spectrum we derived a 
power-law spectral slope /3 — 1.06^013, a soft X-ray absorp- 



+ 13.9 



x 10 2 



Chandra et all J2008l ). IUpdike et al.l J2008I ) and 
20ld ). covering radio to optical frequencies. We 



tion at the host galaxy redshift, A r (H) = 15.5 
cm~ 2 and xt = 0.7 for 40 degrees of freedom (dof). The 
quoted errors correspond to 90% confidence. The mod- 
elled absorbed flux (0.3-10 keV) is F = 1.83 x 10~ 13 ergs 
cm -2 s _1 . The uncertainty on the host N(H) is sufficiently 
large that it would be consistent with only Galactic fore- 
ground absorption. Given that the Galactic iV(H) is typi- 
cally 10 2 ° < N{R) < 10 21 cm" 2 (with uncertainties of ab out 
~ 10% (|Dickev fc Lockmanlll990l ; iKalberla et"aLll2005l ), it 
is often not possible to tightly constrain a low GRB host 
galaxy JV(H). 



3.2 Optical spectrum 

The absorption lines detected in the FORS spectra are listed 
in Table [2] along with the equivalent widths (EWs) . The 
absorption lines associated with the z = 1.5477 system are 
shown in Fig. [2] The apertures used to measure the EWs 
are derived from the resolution in each grism, i.e. twice the 
full width at half maximum of the arc lines (2 x FWHM^rc), 



3 The three GCN measurements do not strongly influence the 
NIR light c urves, since they are v ery close to the NIR photometry 
reported bv lUpdike et al.l 1 1200a) at 2.9 after the burst. 
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Table 1. VLT/FORS2 observation log on date 2007 January 26. a Mid-exposure time after the burst trigger (07:20:45 UT). 



Start time 


At a 


Exp. time 


Grism 


Coverage 


Resolution 


Airmass 


Seeing 


S/N /pixel 


UT (hh:mm:ss) 


(hrs) 


(min) 




A (A) 


FWHM 


(A) 




tt 


(mean) 


04:25:24 


21.20 


15 


300V 


3190-9640 


10.7 




1.8 


0.80 


54 


04:53:31 


21.80 


30 


600z 


7460-10740 


5.3 




1.8 


0.72 


39 


05:25:14 


22.32 


30 


1400V 


4630-5930 


2.0 




1.9 


0.87 


39 


05:56:40 


22.85 


30 


1200R 


5870-7370 


2.2 




2.1 


0.62 


43 



justified by the narrowness of both the Mgn and C IV lines. 
These lines are unresolved, the width of the Mg II A 2803 A 
line being FWHM 2 $ 3 ~ 70 km s _1 . 

The low spectral resolution does not allow us to di- 
rectly derive the column density from the detected EW in 
the optically thin limit, even for the weak EWs we mea- 
sure, since the lines may be saturated. In particular, this 
could affect the Mgn, despite the EW2796 / EW2803 ~ 2 ra- 
tio of the Mgn A 2796 and Mgn A 2803 A suggests little or 
no sa turation. A proper curve of growth (CoG, e.g. jMihalas! 
Il978f ) analysis would require many transitions to constrain 
the 6 value (Doppler parameter, b = \/2 x (TGaussian). How- 
ever, we can benefit from the excellent S/N and the nar- 
rowness of the transitions to constrain the column den- 
sity with the CoG, despite the limited number of detec- 
tions. We studied the CoG evolution with the b value of 
the Mgn AA 2796, 2803 transitions. The rest-frame EW of 
the absorption lines fall between the linear and the satu- 
rated part of the CoG, i.e., the Mgn column density has a 
slight dependency on the b value. We can limit the plausi- 
ble range to 10 < b < 44 km s _1 , as the line broadening 
must be narrower than the observed FWHM ~ 70 km s _1 
(FWHM = 272 In 2 b/^2 = 1.665 b). The lower limit is se- 
lected from the lowest velocity width of low-i onization line 
profile s observed in UVES QSO absorberQ i|Ledoux et al.l 
l2006a ). The A 2803 absorption line gives the tightest con- 
straint, log AT (Mg 11) = 12.96155:18 ( lcr errors), derived con- 
sidering the median column density within the b and lcr EW 
interval. The Mg 11 A 2803 CoG is shown in Fig. El So far 
we used EW2SQ3 and the limits on b to derive the Mg II col- 
umn density. On the other hand, the Mgn column density 
derived from the A 2803 allows us to use the A 2796 CoG to 
further constrain the b value to be 30 < b < 44 km s _1 . With 
these limits on b we cannot derive the C IV column density, 
because the A 1549 CoG shows that the absorption line is 
saturated. We list the 3cr upper limits on the non-detected 



4 The lower observed value is Av = 17 km s —1 , where Av is the 
velocity width that covers 95% of the complex line profile. If the 
line profile is a single Gaussian then Av = 1.96 <7G auss i an = 1.39 b, 
corresponding to a minimum observed b = 12 km s — x . This rep- 
resent only an estimate of the width of the observed line pro- 
file, despite the single components of the complex profile not be- 
ing resolved with low resolution spectroscopy. The UVES spec- 
trum of GRB 050730 showed single component line profiles with 
b = 10 km s _1 . Of course, the possibility that the line pro- 
file is a single component with a very low b can never be to- 
tally excluded. This would result in higher saturation of such 
component, leading to a higher value of 7V(Mgn). A higher 
A(Mgll) would only strengthen the conclusions of this paper. 
The Mgn column density derived in the optically thin limit is 
logJV(Mgll) = 12.88 j 
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Figure 2. Absorption lines detected in the FORS spectra at host 
galaxy redshift 2 = 1.5477. Note that the Civ spectrum, ob- 
served with the lowest resolution, spans a larger velocity range 
(top panel). The error bars show the mean error spectrum, av- 
eraged within ±2000 km s -1 from the line position. The telluric 
features are highlighted. 



Table 2. Detected absorption lines in the VLT/FORS spectra. 
a la errors. Significance 



Ion 


Arest (A) 


-W ( A ) 2 


EW? cst (A) 


a b 


300V 


Civ 


1549 


3949.55 1.5497 


0.77 ±0.09 


8.5 


1200R 


Mgn 
Mgn 


2796.35 
2803.53 


7124.30 1.5477 
7142.74 1.5478 


0.29 ±0.04 
0.14 ±0.02 


7.2 
7.0 


1400V 


Nal 
Nal 


5891.58 
5897.56 


5891.58 0.000 
5897.56 0.000 


0.68±0.08 

0.41±0.08 


8.5 

5.1 



3+0.06 
-0.07' 



lines in Table [3] measuring the noise level at each line posi- 
tion, again with 2 x FWHM^ TC apertures. The excellent S/N 
allows us to put strict constraints on the column densities 
of some ions, such as Fen. 

We identified Galactic (z — 0) absorption of the Nal A A 
5891,5897 doublet in both 1200R and 1400V FORS spectra 
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Figure 3. Curve of growth for Mgn A 2803 for different b values. 
The dotted lines correspond to the measured EW TC st- Including 
the la errors on the EW, we derive logiV(Mgn) = 12.96±o jg. 



Table 3. 3a line detection limits with corresponding column den- 
sities computed in the optically thin limit. 



Ion 


Arost (A) 


EW vcst (A) 


log N 


Grism 


Cn 


1334.53 


< 2.85 


< 15.15 


300V 


Si IV 


1393.75 


< 1.91 


< 14.32 


300V 


Sin 


1526.71 


< 0.44 


< 14.23 


300V 


Alii 


1670.79 


< 0.29 


< 12.80 


300V 


Nin 


1741.55 


< 0.27 


< 14.37 


300V 


Si 


1807.31 


< 0.25 


< 13.89 


300V 


Aim 


1854.72 


< 0.09 


< 12.76 


1400V 


Tin 


1910.93 


< 0.07 


< 13.36 


1400V 


Znll 


2026.14 


< 0.06 


< 12.55 


1400V 


Cm 


2056.25 


< 0.07 


< 13.23 


1400V 


Fen 


2382.76 


< 0.07 


< 12.61 


1200R 


Fen* 


2396.35 


<0.07 


<12.64 


1200R 



and we report the EW in Table[2] The Galactic absorption is 
shown in Fig. IA1I The Galactic reddening in the direction o f 
the target measured by Schlegel, Finkb einer fc Davis! (| 19981 ) 
is E(B — V) — 0.052 ±0.008 mag. For comparison, we derive 
the redd ening from the Galactic N a I - E(B — V) correlation 
given bv lMunari fc Zwitterl IJ1997T ). finding a higher redden- 
ing value, E(B — V) ~ 0.20 mag, but consistent within the 
reddening uncertainty in their correlation (0.15 mag). 



derived by interpolating or extrapolating the light curves 
where necessary, for each filter. The flux interpolations at 
particular times were performed locally, i.e., using up to 
three adjacent data points on each side. The X-ray spectra 
at the five epochs were produced by scaling the overall 
X-ray spectrum to the actual count-rate at each time of 
interest. This choice allowed to preserve the best S/N, 
justified by the lack of spectral evolution measured for the 
X-ray data alone (see Sec. l3.ip . 

The SED was created and fitted in count space using the 
ISIS spectral fitting package (Houck & Dcnicola 2000) and 
following the method given by IStarling et al.l (|2007l ). hav- 
ing the advantage that no model for the X-ray data needed 
to be assumed a priori. Furthermore, as all the data are 
flux-calibrated, no additional offset were needed. We fit the 
SED using models consisting of an absorbed power law or 
absorbed broken power law with slopes free or fixed to a 
slope difference of A/3 = 0.5, as expected for a cooling break 
l|Sari. Piran fc Naravanlll998i ). The latter model is referred 
to as "tied broken power law", since the spectral slope at 
lower energies than the break energy E^, /3i, is tied to 
the spectral slope at higher energies, fa (/3i = /?2 — 0.5). 
The intrinsic optical extinction was modelled with either 
Milky Way (MW), Large Magellanic Cloud (LMC) or Small 
M agell a nic C loud (SMC) extinction curves, as parametrized 
by IPeil |l992j) . The X-ray absorption is modelled with the 
XSPEC model zphabs, assuming Solar metallicity. We also 
tested the model for a lower metallicity. The results for 
each independent epoch are presented in the appendix, Ta- 
ble IA2I We n ote th at the flaring activity, also reported by 
I Up dike et al.l (120081 ) . significantly affects the SED modelling 
at early epochs. 

We also explored the SED evolution by tying some pa- 
rameters between the epochs, where they are not expected to 
change in time. The tied parameters are the host E(B — V) 
and A r (H), since we do not expect the host extinction and 
absorption to vary at two days after the burst, and the spec- 
tral slope at higher energies, /?2, since we detected no signif- 
icant X-ray spectral evolution ( Sec. 13.1]) . The spectral slope 
at lower energies, /3i, and E^k were left free to vary between 
the epochs. We excluded epoch I and epoch II from the SED 
fit tied between the epochs to avoid the flaring activity at 
early times, see Fig. [T] These joint fits are listed in Table 0] 
and Fig. [4] shows the best fitting joint model to the last 3 
epochs. We discuss the best fitting model and a comparison 
with the results of broadband SED modelling in Sec. \S\ 



4 AFTERGLOW MODELLING: THE 
NIR-TO-X-RAY SED 

In order to investigate any temporal evolution of the 
spectral energy distribution (SED) from the NIR to X-ray 
frequencies, we evaluated the SED at five epochs after the 
burst trigger. The SED times 1.0, 1.4, 1.9, 2.5, 2.9 days 
were selected in order to maximize the NIR to UV data 
coverage for each epoch and minimize extrapolations, as 
shown in Fig. [l] For each SED epoch, we combine the 
available NIR to UV photometric data points and the ones 



5 DISCUSSION 

5.1 The spectral energy distribution 

The jointly fit epoch III-V NIR-to-X-ray SED of 
GRB 070125 is best described by a broken power law with 
/?2 = 1.18 ± 0.01 (tied between the epochs), and /3i = 
0.835±g;g||, Q-72±g;g|, 0.45±°;^ for epoch III, IV and V re- 
spectively (xt = 2.84 for 146 degrees of freedom, dof, see 
Table [4} . The break energy, i?bk is located in the U V and 
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Figure 4. Broken power-law fit to the epoch III (blue), IV (cyan) and V (magenta) NIR-to- X-ray SED, with LMC dust extinction and 
tied parameters between the epochs. 



there is a marginal hint for its evolutioqj in time towards 
longer wavelengths: A bk = 1632±|gl, 2340^, 2885^ A 
from epoch III, IV and V. The quoted errors refer to 90% 
confidence. The spectral slope differences suggest that the 
spectral break is a synchrotron cooling break, and a fit with 
a tied broken power law (A/3 = 0.5) provides an adequate fit 
indeed (\l = 2.97 for 147 dof) and the same /3 2 = 1.18±0.01. 
The free broken power law is slightly preferred with respect 
to the tied broken power law, but the latter cannot formally 
be excluded (F-test null probability of 0.62). Regardless of 
the extinction model used, a broken power law always gives a 
significantly better fit than a single power law from the NIR 
to X-ray regime. The F-test probability of the broken power- 
law over the power-law fit improvement being obtained by 
chance is between 10 -12 and 7 x 10~ n with all extinction 
laws in the epochs III- V joint fits, confirming the significance 
of the improvement. 

The dust extinction A v - 0.05 ± 02 m ag (for R v - 
A V /E(B - V) = 3.16 in the LMC, |Pelll992l ) is best mod- 
elled with a LMC or SMC extinction curve, suggesting a 
low metallicity GRB absorbeiQ We note that the free bro- 
ken power law with MW-type extinction curve cannot be 



5 Consistent with a time evolution of the synchrotron cooling 
frequency as u c oc t~ 1 ' 2 expected for a blast wave expanding 
in a co nstant ISM medium in the s low cooling regime [y m < 
u c , e.g., ISari, Piran fc Naravanl|l998l ). However, the evolution is 
not significant, since fitting a constant to the three E^k values 
provides x 2 = 1-4 (for 2 dof). 

6 [Fe /Hl = -1.2 5 ± 0.01 ( 1.15) in SMC (LMC) star clus- 
ters llCionil 120091) . while the Milky Way metallicit y distribution 
peaks around Solar values both in the bulge (e.g., IZoccali et alj 
2003) and in the thin disk (Tiede & Terndrup 1999), covering the 
majority of the stellar population. 



formally excluded by the fit (xl, — 2.99 for 146 dof). How- 
ever, the very low dust content implied by this solution 
(E(B - V) = 0.008 ± 0.007) also suggests a low metallic- 
ity for the absorbing gas (e.g.. lVladild ll998). 

7V(H) cannot be constrained, i.e., it is consistent with 
the Galactic value and with the large error of the X-ray- 



only JV(H) = 15.5 



+ 13.9 



x W 



This JV(H) falls at 



the l ow end of the -W(H) distribution of GRB afterglows 
(e.g., lEvans et al.|[2009l ; iFvnbo et al.ll2009l ; ICampana et all 
[2010J), and Av also falls at the low end of the dust ex- 
tinct i on distribution ( e.g., iKann et al.l 120101 ; ISchadv et al.l 
l201ll ; IZafar et alj|201lt ). The dust to gas ratio N(K)/A V = 



3.1 



+2.7 



x 10 cm mag covers most of the GRB 



N(H.)/Av distribution at its lo w end (seelZafar et al.ll201l|j 
and is similar to GRB 980329 ^Starling et al.ll2007l ). In fact, 
one in every six Swi ft XRT X-ray afte rglows have uncon- 
strained excess iV(H) (|Evans et alj|2009f ). given the difficulty 
of measuring low host galaxy iV(H) against higher Galactic 
iV(H) and its uncertainties. 

At earlier epochs I and II the afterglow SED is affected 
by flaring activity. This is most evident in the V and R bands 
where epochs I and II coincide with the rise and decay of a 
flare or bump in the light curves (Fig. [TJ . For completeness 
we report individual fits to these early epochs in the ap- 
pendix. Our NIR-to- X-ray SED results are in overall agree- 
ment wi th the single epoch S ED at 4.3 days previously anal- 
ysed by lUpdike et al.l (|2008l ). In particular, we confirm the 
presence of a spectral b r eak. T his break was only tentatively 
found by lUpdike et alj 1)20081 ). because of its sensitivity to 
the NIR continuum level, based on a large extrapolation of 
solely two data-points in the K band. In our dataset, 6 ad- 
ditional NIR data-po ints were added, i .e. K s , J,H from the 
SMARTS telescope (|Kann et all l2010h and K S ,J,H from 
the PAIRITEL telescope (from GCN. IChandra et al-lliooih . 
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Table 4. NIR-to-X-ray SED modelling with a single power law (PL), a broken power law (BPL) and a tied br oken power law (TIED 
BPL) , where the spectral slopes are tied to differ by A/3 = 0.5, as expected for a synchrotron cooling break iSari, Piran fc Naravanl 
1998). Each model is jointly fit to the SED at epoch III, IV and V (1.9, 2.5 and 2.9 days after the trigger) by tying some parameters 
between the epochs, i.e. E(B — V), N(H) and 02 are allowed to vary but not between the epochs. The best fitting model (underlined) 
is a broken power law with 02 = 1.18 ± 0.01 and a break energy E h 



7. 



3+4.0 



5.3 



+ 1.2 
1.0' 



4.3+2° eV ' in the UV ' The iV ( H ) is consistent 
with the Galactic value, while a low dust content E(B — V) = 0.016 ± 0.007 mag can be derived from the SED, equally well described 
by SMC and LMC extinction curves. The errors refer to a 90% confidence limit, and are solely statistical error, i.e. they do not include 
systematic errors. 



Ext. 
type 


Model 


Xl [dof] 


E(B-V) 
(mag) 


iV(H) 

(10 22 cm" 2 ) 


01 


02 


(eV) 


Epoch 




PL 
BPL 
TIED BPL 


4.06 [148] 
2.85 [146] 
2.96 [147] 


< 0.0018 
0.015 ±0.007 
0.017 ±0.006 


< 1.3 x lO" 6 

< 1.2 x 10~ 6 

< 1.2 x 10~ 6 


1.083 ±0.008 
1.091 ± 0.009 

1 069+ 0010 

l.UD9_ 009 

0.828 ±0.065 

n 7 9+0.08 
u -' -0.09 

0.48± , : 1 2 3 

02 - 0.5 






III 
IV 
V 

III 
IV 
V 

III 

IV 
V 














SMC 


1.18 ±0.01 
1.18 ±0.01 


7 5+3.7 
'•°-1.8 
5 4+1-3 
°- 4 -0.9 

4 5+0-9 
4 -°-0.7 

5 1+0.9 

51+O.8 
°- 1 -0.7 

6 0+ 1 - 1 
°- u -0.8 




PL 
BPL 
TIED BPL 


4.06 [148] 
2.84 [146] 
2.97 [147] 


< 0.0018 
0.016 ±0.007 
0.016 ±0.006 


< 1.3 x 10~ 6 

< 1.3 x 10~ 6 

< 1.2 x 10~ 6 


1.083 ±0.008 
1.091 ± 0.009 

1 069+ - 010 
i.uoy_ 009 

n aqrc+0.065 
U.OOO_ os3 
n 7 9+0.08 
u -' -0.09 
(,45+0.12 

02 - 0.5 






III 

IV 
V 

III 

IV 
V 

III 

IV 
V 














LMC 


1.18 ±0.01 
1.18 ±0.01 


76 +4.0 
'• u -1.8 
5 o+1.2 
o - o -1.0 

4 s+ - 9 

5 9+0.9 
°- z -0.7 

5 1+O.8 
a -l-0.7 

6 0+ 1 - 2 
°- u -0.8 




PL 
BPL 
TIED BPL 


4.06 [148] 
2.99 [146] 
4.14 [147] 


< 0.0017 
0.008 ± 0.007 

< 0.0016 


< 1.3 x 10~ 6 

< 1.2 x 10~ 6 

< 1.2 x 10~ 6 


1.083 ±0.008 
1.091 ±0.009 
1 069+ 0010 

l.UD9_ 009 

1 n8Q+ 0010 
l.uoy_ 009 

o.ni- ;°' 
0-50 + °;i° 

02 - 0.5 






III 

IV 
V 

III 

IV 
V 

III 

IV 
V 














MW 


1.184 ±0.016 

i rcoi +0.006 
l.OOl_ 005 


6 1+ 1 ' 7 
°- 1 -i.o 

< 5 

10.00lg;°g x 10 3 
3.17to'?7 x 10 3 
7.58t 2 j 2 x 10 3 



Moreover, no extrapolations were needed for the NIR light- 
curve from epoch III to epoch V of our SED fit. Finally, 
the significant improvement from a single power law SED 
fit to a broken power law strongly favours the presence of 
the spectral break. 



5.2 The broadband picture 

To determine the full set of physical parameters of a GRB 
afterglow, it is vital to have good temporal and spectral cov- 
erage from X-ray to radio frequencies. GRB 070125 is a good 
candidate for this kind of study, in particular because of its 
coverage at centimetre wavelengths. With well sampled light 
curves at multiple radio observing frequencies one can de- 
termine the evolution of the peak flux, the associated peak 
frequency i^ m , and the synchrotron self-absorption frequency 
v & . Together with the cooling frequency u c , the electron en- 
ergy distribution index p, and the jet-break time, all the 
micro- and macrophysical parameters of the GRB jet and 



its surroundings can be pinned down. F or GRB 070125 this 
has been done by I Chandra et al.l (120081) for part of the total 



available data set. I Up dike et al.l IJ2008I ) have focused their 



modelling efforts mainly on the optical and X-ray regimes, 
as we have done in Sec. 15.11 There are some inconsistencies 
between the different modelling results, and complications 
in interpreting the full data set, which we will discuss here. 

The main issue is the position of the cooling break in 
the broadband spectrum. Our SED fits show that there is 
a spectral break between the optical and X-ray bands. The 
spectral slope differences between the two bands indicate 
that this break can be interpreted as the cooling break. As 
pointed out bv lUpdike et al.l (|2008f ). it is difficult to form a 
coherent picture of these spectral slopes and the temporal 
slopes observed at optical and X-ray frequencies within the 
standard afterglow fr amework, because not all the so-called 
closu re relations (e . g ., Zhang fc MeszarosH2004T ) can be satis- 
fied. ly_piiike_et_alj (J2008I ) suggested a solution for this issue 
by adopting that the NIR-to-X-ray SED is in fact a sin- 
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gle power law, and they concluded that all those observing 
bands are in between z/ m and u c and that the circurnburst 
medium is homogeneous. This solution is not satisfactory 
given that we confirm the presence of a cooling break in the 
NIR-to-X-ray SEP , similar to what w as found and tenta- 
tively discarded bv lUpdike et al.l (120081 ). 

A complicating factor is the flaring behaviour of 
GRB 070125, which makes reliable determinations of tem- 
poral slopes difficult. Our SEDs have been constructed at 
times when the flaring has ceased, but leaving out the flaring 
parts of the optical and X-ray light curves gives a relatively 
short lever arm for determining the temporal slopes, in turn 
giving rise to large statistical and systematic uncertainties. 
This could account for the fact that the standard closure re- 
lations for u c in between the optical and X-ray bands do not 
seem to work. For example, it is mentioned in I Up dike et al.l 
(2008) that the optical temporal slopes are consistent with 
the spectral slopes if the circurnburst medium is structured 
like a stellar wind, but the X-ray temporal slope is incon- 
sistent. The latter could be due to uncertainties caused by 
the flaring of the source. The bottom line is that in the case 
of GRB 070125 the NIR-to-X-ray SED fits seem to give the 
most reliable measurement of the value of v c . 

The bro a dband SED fitting performed by 
IChandra et al.l l|2008l ). including the centimetre wave- 
length radio data, leads to different conclusions. They 
find that the optical and X-ray regimes lie on the same 
power-law segment of the synchrotron spectrum, but above 
both Vm and v c . In fact, one of the outcomes of their 
modelling is that this afterglow is in the fast cooling regime, 
i.e. v c < z/ m , for about a week, which has not been observed 
in any GRB afterglow before. The values for p they find 
result in optical and X-ray spectral indices of ~ 1.1 — 1.2. 
These indices are consistent with the X-ray spectral index 
we find in our SED fitting, but inconsistent with the shal- 
low er optical indi c es we find, and also the indices obtained 
by lUpdike et all (|2008l ) for NIR-optical SED fits. The 
uncharacteristically low value for u c is necessary to describe 
the optical and X-ray light curves together with the bright 
and long-lasting radio afterglow. If v c would be at higher 
frequencies, there is no possible combination of z/ a , v m and 
the peak flux to give a satisfactory fit of the radio data . 
Although the modeling performed bv lChandra et al.l (| 20081 ) 
is thorough, the basic assumptions are fairly simplified. 
To get to a coherent solution for all the broadband data, 
a more complicated picture may have to be invoked, like 
a double jet as has been s uggested for the b roadband 
afterglow of e.g. GRB 030329 l|Berger et al.l 120031 ). Another 
possibility is that we need a more detailed description of the 
afterglow physics, in particular the effects of synchrotron 
self-absorption on the jet-break time at radio wavelengths 
and the associated chan ges this makes to radio light curves 
|van Eerten et al.ll201ll ). These different, more complicated 
scenarios for satisfactorily describing the full data set of 
GRB 070125 are beyond the scope of this work, but will be 
explored and presented in a future paper. 

While the broadband modelling needs further study, we 
can still rely on the NIR-to-X-ray SED, as it is based on 
fewer assumptions. In particular, we found a low iV(H) and 
a low SMC/LMC-type dust content, suggesting a low metal- 
licity environment. 



5.3 The absorber environment 

While the broadband SED depends on the afterglow physics 
and the circurnburst environment, the absorption in the op- 
tical spectrum probes the gas along the line of sight. By 
studying the nature of the absorber we can derive informa- 
tion on the GRB site, as we will discuss here. 

A crucial issue that needs to be addressed is whether 
or not the Mg II absorbing system observed in the afterglow 
spectrum is associated with the host galaxy of GRB 070125. 
While a host galaxy could, in principle, be in the background 
at a larger redshift, we note that Lya absorption at z > 1.9J 
would be observable in our high S/N FORS spectra, but no 
signature of such system was detected. Figure [5] shows how 
the observed GRB 070125 FORS spect rum compa r es to the 
GRB composite spectrum from lChristensen et al.l (|2011f ) at 
the same redshift z — 1.5477. A 1.5 < z < 1.9 host galaxy 
could be possible but no lines are detected at such redshift. 
Thus, we consider z — 1.5477 to be the host galaxy redshift. 
A tentative Lya abso rption at this redshift was reported 
bv lUpdike et al.l (|2008l ). for which log AT (Hi) < 20.3 was ar- 
gued. We note that high-ionization N v AA 1238, 1242 A lines 
are not in our wavelength reach. No spectroscopic observa- 
tion has been attempted to inves tigate the two candidate 
blue hosts that ICenko et al.l (|2008i) indicated at a large off- 
set from the GRB position so far. Below we will discuss the 
possible scenarios. 

The metallicity along the GRB 070125 line of sight can- 
not simply be constrained from the measured column densi- 
ties or upper limits, since no Lya is available in our FORS 
spectra, given the redshift of this afterglow. One way to 
bypass the lack of iV(Hl) measurement is plotting [X/Hl] 
+ logJV(Hi) = logiV(X) - log7V(X) + logiV(Hi) , as in 
Fig. HJ] Despite the degeneracy between Hi column den- 
sity and metallicity, we can use the Mgn detection and the 
upper limits to explore the parameter space and at least 
limit the possible iV(Hl) and [X/Hl]. In particular, assum- 
ing Solar metallicity and no Mg depletion into dust, we 
derive logiV(Hl) = 17.43. If the metallicity is 1/100 So- 
lar then logiV(Hl) = 19.43. A much lower metallicity is 
very hard to justify, while a super-Solar metallicity with 
lower iV(Hl) is possible. Hence, we can reduce the possi- 
bilities for the GRB 070125 absorber to be either i) a sub- 
Lyma n-Limit System (sub-LLS, logiV(Hl) < 17.3. iRauchl 
I1998T ) with Z > 1.3 Z G) , o r ii) an LLS (17.3 < logA^Hl) 
< 19.0. iPeroux et alj|2003l ) with 0.03 Z e < Z < 1.3 Z e , or 
Hi) a sub-DLA (19.0 < logiV(Hi) < 20.3) with 0.001 Z e < 
Z < 0.03 Zq. Below we discuss these three different scenar- 
ios. 

i) A super-Solar Z sub-LLS. This scenario is supported 
by the low Mgn EW of GRB 070125, since the majority 
of weak Mgn (EW < 0.3 A) absorbers in quasar line of 
sights are associated with su b-LLS with super-Solar metal- 
licity l|Churchill et al.lll999l ). This association is based on 
the number density consistency between the two popula- 
tions and a CLOUDY photo-ionization model that predicts 
log N(Hl) < 17.3 for weak Mgn. Given the observed Mgn 
column density in GRB 070125, the C II, Si ii and Fe ii expec- 
tation for weak Mg II systems derived by iNaravanan et al.l 



7 Corresponding to the bluest observable wavelength A = 3515 - 
in the spectrum cut where S/N~ 10 
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Figure 5. The high S/N GRB composite spectrum (dark) from Christcns en et alj (I201lh shifted to z = 1.5477 and the observed 
GRB 070125 FORS spectrum (light) cut at S/N > 10. The error spectra are shown by the dotted lines. A firm upper limit of z < 1.9 can 
be placed on the GRB 070125 redshift from the lack of Lya forest lines in the afterglow spectrum. 
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Figure 6. Abundances and 3<r detection limits of a selection 
of low ionization species in the GRB 070125 afterglow spectrum 
(black bar and arrows) from Table [2] and Table [3] The ordinate 
combines the metallicity and the Hi column density, basically a 
measure of the metal abundances with respect to Solar, in or- 
der to show the possible ranges of metallicity and ./V(Hl) in the 
absorber. The horizontal dashed lines correspond to the neutral 
hydrogen column density for a given metallicity, log N(H l)z, with 
respect to the Mgll detection, assuming no depletion of Mg into 
dust. Light- coloured square s and show the meta l abundances ob- 
served in GRB absorbers (S chadv et al,ll201lT) and the crosses 
refer to the GRB 080310 measurements (De Cia et al., in prep.). 
The column density expectations (diamonds) for the measured 
Mg II ar e derived from the corr elations observed in weak Mg ab- 
sorbers JNaravanan et al.l l2008). 



(2008) (diamonds in Fig. [6]), are consistent with our up- 
per limits. These systems could trace either interstellar gas 
expelled from star forming regions into the halo, or small 
clouds in dwarf galaxies or intergalactic star-forming struc- 
tures (|Churchill et al.lll999l ). The l east likely po s sibilit y is 
the GRB explosion site in the halo. ICenko et al.l (|2008T ) al- 
ready pointed out the difficulty of explaining how to "kick" 
the massive GRB progenitor star all the way out to the 
galaxy halo (~ 27 kpc). A halo in-situ formation of the 
GRB would require halo gas at ~ 27 kpc and ~ 46 kpc 
from two interacting galaxies to trigger massive star for- 
mation, the GRB progenitor producing high density stel- 
lar wi nds to "fuel" the extr emely bright afterglow (n ~ 50 
cm" 3 . IChandra et al.l 120081 ). One oth er possibility i s that 
the GRB exploded in a dwarf galaxy. Ell ison et al.l ()2004l ) 
studied the Mgll systems coherence length, indicating that 
the weak Mgll absorbers have sizes of 1.5-4.4 kpc and they 
could represent a distinct population of smaller galaxies. 
However, the low content of SMC/LMC-type dust derived 
from the GRB 070125 afterglow SED suggests a low metal- 
licity of the host galaxy, disfavouring the super-Solar Z sub- 
LLS scenario. 

n) A 0.03 Z Q < Z < 1.3 Z Q LLS. Although GRB ab- 
sorbers show, on aver age, high Hi column densities (e.g., 
Ijakobsson et al.l 120061 ). an LLS is not an exotic environ- 
ment for a GRB, since a low N(Hl) may be a line of 
sight effect. In principle, the burst radiation can travel 
along the line of sight, within the host galaxy or from 
the galaxy outskirts, without being absorbed by dense gas 
clouds. This scenario is supported by the observation of 
a few other GRB LLSs, 
= 17.60 ± 



f GR B LLSs, namely GRB 050908 (log AT (H i) 
0.10: iFueazza et aD 120051 : iFvnbo et all 1200 9) . 



Prochaska et al.l 



GRB 060124 (logJV(Hl ) = 18.5 ± 0.5: 
| 2006l:lFvnboet al.ll2009t), GRB060607A [log AT (Hi) ~ 16.8: 
Ledoux et all l2006bl : IProchaska et al.l l2008r i. GRB 080310 
(logAT(Hi) = 18.70±0.10: De Cia et al., in preparation) an d 
GRB090426 (logAf(Hi) = 18.70+ijo: pThone et all l20li1 ). 
where Lya variability was induced by Hi ionization. The 
low H I content of all these systems cannot be explained by 
ionization alone, as the burst radiation can totally ionize 
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H I only out to 200-300 pqj from the burst i n LLSs (using 
the b right GRB 050730 afterglow light curve, iLedoux et al.l 
2009). A low Hi could arise if the GRB is not too deeply 
embedded in its host galaxy. Such a possibility, which is cer- 
tainly expected in some cases, will be further investigated 
later on in this Section. 

The equivalent hydrogen column density measured from 
the soft X-ray absorption can also provide important infor- 
mation on the gas along the line of sight, with the advantage 
of not depending on th e Hi ionization. Indeed, iV(Hl) and 
JV(H ) do not correlate l|Watson et al.ll2007l ; ICampana et alj 
120101 ) . The difference between a low neutral and a high total 
hydrogen column densities can be explained if most of the 
X-ray absorbing gas lies closer to th e burst, where the gas 
is highly ionized i Schadv et al.ll201ll ). In this scenario, the 
GRB 070125 could still have been surrounded by the high 
density gas required to power such a bright afterglow. 

Hi) A 0.001 Z Q < Z < 0.03 Z Q sub-DLA. This metal- 
licity range is low, but not unusual for GRB absorbers. For 
instance, most of the GRB absorbers in the VLT/UVES 
sample have 0.006 Zp, < Z < 0.05 Z q, despite the small 
number statistics (|Ledoux et al.l 12009 ). although their Hi 
column densities are generally higher. 

From the above analysis, scenarios ii) and partly Hi) 
are favoured, i.e., the GRB 070125 environment is consistent 
with other GRBs. Below we test this further, by comparing 
the GRB 070125 properties with the GRB population. 

In Fig. [6] we compare the measured metal column den- 
sities and upper limits for GRB 070125 with those observed 
in other GRB aft e rglow s (squares and arrows), as collected 
by ISchadv et al.l (|2011f ). The metal column densities in 
GRB 070125 are generally lower than measured in other 
GRBs. This suggests that the GRB 070125 absorber has ei- 
ther a lower metallicity or a lower H I column density with 
respect to most GRBs. However, this comparison sample 
is limited, since the column densities often cannot be con- 
strained, due to the low resolution of the spectra, and more 
complete EW distributions should be used. The EW of C IV 
A 1549 in GRB 070125 is EW rcst = 0.77±0.09 A, much lower 
t han in the GRB composite spectrum (EW TCS t — 2. 18 ±0.03 
A lChristensen et al.ll201ll) and in the low H i afterglow sub- 
sample (EWrest ~ 2.8 A. lThone et all 1201 ll ~). Nevertheless, 
the observed Civ A 1549 EW and the Sin A 1526 up- 
per limit are consi stent with the EW correlation found in 
iFvnbo et ail (|2009l V This suggests that the GRB 070125 ab- 
sorber fits in the GRB population, rather than being an 
outlier. In Fig. [7] we compare the Civ AA 1548, 1550, Mgii 
A 2796 EWs and the deep Fen A 2382 3<r detection limit 
for GRB 070125 with the EW distributions discussed by de 
Ugarte Postigo et al. (in preparation) fo r the whole low reso - 
lution spectroscopy afterglow sample of lFvnbo et al.l (|2009h - 
The GRB 070125 absorber falls in the low end of the EW 
distribution. Again, this indicates that GRB 070125 is not a 
GRB population outlier and that other bursts share similar 
properties. 

Furthermore, we selected two afterglow sub-samples 
that may potentially be similar to GRB 070125: a) the 
spectral-line poor, i.e., afterglows with no more than two 



8 For much brighter afterglows, these distances can significantly 
grow, see our discussion below. 




Figure 7. The Civ A 1548/1550, Mgn A 2796 and Fen A 2382 
restframe EW distributions for GRB absorbers. The EW and 3cr 
detection limi ts are taken from de Ugarte Postigo et al. (in prep.) 
and based on IFvnbo et alj feOOSft . The highlighted GRB 070125 
EWs fall at the low end of the EW distribution. 



species in absorpti on at the GRB redshift, selected from 
IFvnbo et alj i|2009l ) and b) afterglows with low equivalent 
hydrogen column density iV(H) measured from the X-rays, 
N(K) < 5 x 10 21 cm" 2 withi n 3a. JV(H) and A v were col- 
lected from Zafar ct al. (2011) when possible, o therwise the 
N(H) from lEvans et al.l (J2009) and the Ay fro mlKann et al 



(2010) were used. The iV(Hl) were taken from lFvnbo et al 



(2009) and the EWs from de Ugarte Postigo et al. (in prep, 
with the exception of GRB 060708 and GRB 080310, whose 
metal EWs are measured here. We exclude the low S/N 
spectra from the spectral-line poor sample by applying a cut 
in the Mgn A 2796 detection limit of uew < 0.5. The low 
S/N spectra are also excluded from the EW constraints of 
the low iV(H) sub-sample. Six afterglows fall in our spectral- 
line poor sub-sample, not including GRB 070125, listed in 
Table [5] while fifteen are in the low iV(H) sub-sample. 

We note that the spectral-line poor sub-sample may suf- 
fer from selection effects. First, afterglows with low redshift 
[z ^ 0.5) may more easily be included in this sub-sample be- 
cause of the shortage of observable transitions (other than 
the Mgn AA 2796, 2803 doublet) in the spectral region cov- 
ered by typical optical spectrographs. However, only two 
GRBs have such low z in Table [5] Second, the S/N require- 
ment may bias the sub-sample towards brighter bursts. This 
is irrelevant when comparing the sub-sample properties with 
GRB 070125, the latter being very bright itself, but should 
be kept in mind when considering the spectral-line poor sub- 
sample as a population. A more extended and complete sam- 
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Tabl e 5. Spectral-line poor sample (first group, based on the lFvnbo et alj2009 . sample) and rela ted cases (second group, not included in 
the lFvnbo et alj|2003 . sample). a Afterglow brightness a t 1 day after the bu rst, shifted to z = 1 IIKann et al.ll2010l) . b R, i or z afterglow 
brightness from the acquisition i mages of the spectra I Fynbo et al.| l200al . c No observations av ailab le. d Possible contribution from 
th e SN or the afterg l ow its elf. [ 11 ISavaglio. Glazebrook fc Le Borend d2009tl ; [2 1 IChen et alj [|2009t); [31 lLaskar, Bereer fc Charvl J201ll) ; 
[41 iDella Valle et all feoOrj); [5] iGal-Yam et alj J2006I) ; [6] ICovino et alj J2010I) ; [7] IFynbo et al.l J2OO9 ); [8] Cenko et al. (in prep.); [9] 
IPerlev et all l|200St ): [10] ISparre et all l|201ll ). 
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pie with deeper S/N limits is needed to further analyse and 
possibly exclude these selection effects. 

Fig.[8]shows a comparison of the GRB 070125 afterglow 
absorber properties with the two sub-samples. GRB 070125 
traces the low end of the iV(H), Av and E W distributions, 
confirming the previous finding for the whole sample. 

In general, we notice that most of the spectral-line poor 
afterglows also show a low iV(H). The few low A/(H) after- 
glows with Lya measurements seem to prefer low Hi envi- 
ronment, although no conclusion can be drawn yet due to 
the extremely poor statistics. The same applies to the lack of 
high EW at constrained low iV(H). The above suggestions, 
if confirmed, may indicate a overall low column density en- 
vironment for both the low-ionization gas (that absorbs the 
metal lines in the optical spectrum at distances larger than 
a few hundred pc from the burst) and the highly ionized gas 
(that absorbs the soft X-rays, in the vicinity of the burst). 
More observations are needed to investigate whether N(H) 
and iV(H 1) could correlate under special conditions, e.g. due 
to the effect of a short line of sight within the host galaxy, 
despite there generally being no such correlation. However, 
we note that the A/(H) typically cannot be constrained be- 
low log A/(H) < 20-21, due to the uncertainty in the Galactic 
N(S). 

One GRB absorber potentially similar to GRB 070125 
is the recent nearby GRB 101219B, with redshift 2 = 0.55 
identified from the Mgn AA 2796, 2803 A and Mgi A 2853, 
the only features in the afterglow spectrum. A spectroscop- 
ically confirmed supernova was associated with this GRB, 
supporting both the low redshift and a massive star progen- 
itor for the GRB. The afterglow spectral shape suggested 
a very low amount of dust (Av < 0.1), while the non- 
detection of the host galaxy continuum or emission lines 
indicated that a small galaxy was hosting the GRB , at the 
faint end of the Savaglio. Glazebrook fc Le Borgnd (2009) 
distribution (jSparre et al.ll201ll ). 

Another case analogous to GRB 070125 was 
GRB 071003, a very bright burst (R « 12 mag about 
a minute after trigger), with few and very weak absorption 
lines in the late time spectrum and no underlying bright 
host. Given the similarity between the two, GRB 071003 



was also claimed to have occurred in a halo environment 
IjPerlev et al.ll2008l ). In both cases, while the high brightness 
of the afterglow cannot easily be associated with a halo 
star-fornring region, it can naturally contribute to ionize 
the circumburst medium to larger distances. And perhaps 
the time delay between the burst and the epoch of the 
spectra (~ 1 day in both cases) could allow the ionization 
front of the bright afterglow to travel further through the 
host galaxy. 

In order to investigate this possibility, we performed a 
simplified simulation of the GRB 070125 afterglow radiative 
effect on the surrounding medium with a radiative transfer 
model, including photo-ionization. A detailed description of 
the complete model will be discussed in Vreeswijk et al. (in 
prep). We use the GRB 050730 light curve as a template, 
scaled to the 3.3 times brighter flux of GRB 070125 (derived 
from the afterglow brightness at 1 day, if both GRBs would 
lie at z = 1, .fie (G RB 050730) -Re (GRB 070125) = 1.3 mag, 
iKann et al.ll2010l ). We assume a spectral slope j3 — 0.7, as 
before the synchrotron cooling, since the most relevant pho- 
tons are the UV in the first hours after the burst. We let the 
GRB afterglow radiate through a cloud lying at a distance 
between and 2 kpc and investigate the variation of the 
ionic column densities. In general, most of GRB absorbers 
have been found to lie within this distance range from the 
burst, in those cases for which the GRB-to-cloud distance 
could be estimated from the variability of fine-structure lines 
in the spectra. Two different sets of initial column densities 
were chosen: a) N(Ul) = 10 19 cm" 2 and N = 10 14 cm" 2 
for the metal ions and b) AT (Hi) = 10 20 cm" 2 and N = 10 15 
cm" 2 for the metals. We set the observation times to be at 
1 day and 6 hr for comparison (observer frame). The results 
for case a) are displayed in Fig. [5] While the higher N(H 1) 
shields the metals from being ionized at large distances, in 
the case of N(Ri) = 10 19 cm" 2 , most of the metals are 
strongly ionized out to ~ 0.5 kpc. In particular, at 0.5 kpc 
the Mgn column density is 10 times lower than the initial 
value, Fell is ~30 times lower than the initial column, Nill 
even lower. The late time of the spectrum (1 day) does not 
seem to be an important factor, as the comparison with the 
ionization expected at 6 hr after the burst does not change 
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Figure 8. The neutral hydrogen column density, dust extinction and restframe EWs vs the equivalent hydrogen column density of the 
spectral-line poor (dark) and the low iV(H) (light) GRB sub-samples. The GRB 070125 properties are highlighted. The dotted line in the 
top panel marks where AT(H) and JV(Hl) are equal. We do not display the four highest Ay values (0.34 < Ay < 0.49) at the high iV(H) 
distribution end, in orde r to fo cus on the low Ay and N(R) cases most similar to GRB 070125. A complete iV(H) vs Ay distribution is 
presented in lZafar et al.l |2011f) . 
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Figure 9. The ionization effect of the afterglow radiation on 
absorbing clouds at a given distance with initial column densities 
logA/(Hl) = 19 and logW = 14 for the metals, as observed at 1 
day after the burst. See text for details. 



the result dramatically. In particular Mg n and Fe n are ion- 
ized by only a factor of ~2 between 6 and 24 hours, indicat- 
ing that most of the ionization is happening within the first 
hours after the burst. At 0.2 kpc all these metals and Hi are 
totally ionized. Thus, we conclude that any cloud within 0.5 
kpc would not be visible, in terms of spectral lines of low 
ionization species. This evidence, together with a GRB loca- 
tion not too deep inside the host galaxy (i.e., within 0.5 kpc) 
can explain the lack of spectral features in the spectrum of 
GRB 070125. 

The feasibility of this scenario for other bursts can be 
tested by comparing how many cases similar to GRB 070125 
we expect and how many we actually observe. For a zero or- 
der approximation, we calculate the rough probability of a 
GRB to have i) a short line of sight within the host galaxy, 
ii) a low absorbing iV(Hl) and in) a bright afterglow, i) 
The probability of a GRB to be located not deeper than 0.5 
kpc within its host galaxy depends on the galaxy size (in 
one simple dimension it can be simplified as p(< 0.5 kpc) ss 
0.5kpc/size. Integrating this probability over a normal dis- 
tribution of GRB host galaxy sizes centred around t he mean 
size 2.4 kpc (jWainwright. Berger fe Penprasd 120071 . assum- 
ing a = 0.8 and excluding sizes below 0.5 kpc, to reproduce 
the observed 0.5-5 kpc size distribution) provides a 22% 
chance of a GRB line of sight traversing less than 0.5 kpc 
of its host galaxy. |j ii) Among the GRB absorbers with de- 
tected Lycv. absorption in the iFvnbo et al.l (120091 ) sample, 
15% are LLS or lower AT (Hi). Hi) The afterglow bright- 



ness distribution, scaled to re dshift z = 1 and at 1 day 
after the burst, is presented in iKann et al.l (J2010j), provid- 
ing 16% of GRBs with R c < 18 mag afterglows. A differ- 
ent approach, considering the magnitudes of the acquisition 
images ( R, i or z-ban d mag nitude) of the afterglow spec- 
tra from IFvnbo et al.l l|2009h , yields 12% afterglows with 
Magacq < 18 mag. Combining the three probabilities above 
we estimate that about 0.5% of all GRB afterglows are simi- 
lar to GRB 070125, i.e. spectral-line poor and bright. This is 
somewhat consistent with the observations of GRB 070125, 
GRB 060607A and including GRB 071003, the t hree of them 
repres enting 4% of the afterglows analysed in IKann et al.l 
(|201Ci ) . We made an approximate estimate of the number of 
GRBs with a ground-based detection among the Swift sam- 
ple (and hence the possibility of being identified as spectral- 
line poor, low iV(HI) and optically bright), using the Swift 
GRB Tabl«E3- This results in 324 GRBs up to 2011 May 1, 
of which 1% would mean ~3 sources, again consistent with 
current observations. If we relax the brightness requirement, 
about 3% of all GRB afterglows are expected to have a low 
AT (Hi) and lie in the outer kiloparsec of their host galaxy 
along our line of sight. The 9 spectral- line poor afterglows 
listed in Table [5] (including GRB 71003 and GRB 10 1219B) 
represent 12% of the afterglows of lKann et al.l (|2010j), or 3% 
of the Swift optical/IR afterglows. It is possible that these 
spectral-line poor afterglows are the result of a GRB location 
in the outer regions of their host galaxy and a low 7V(Hi). 
Further investigations are needed to confirm these results. 

Finally, we point out that the non-detection of the 
GRB 070125 host g alaxy in the Keckl /LRIS images down 
to R > 25.4 mag l|Cenko et al.l [2008) is consistent with 
a M300, ab > —18.8 galaxy at z = 1.5477 (k-corrected 
AB magnitude at 300 nm), i.e., a faint (among the 
25% faintest) but still not uncommon GRB host galaxy 
at that redshift (Schulze et al. 2011, in preparation). 
Ind eed, GRB hosts are generally faint ir regular galax- 
ies l|Savaglio. Glazebrook fe Le Borgnell2009h . Recent HST 
imaging of the field provides deeper upper limits on the 
host brightness of F336W > 26.7 mag (AB) (central wave- 
length A c ~ 334.4 nm) and F110W > 26.4 mag (A c ~ 1123.8 
nm), not corrected for foreground extinction (Cenko et al., 
in prep.). We then correct for Galactic foreground extinc- 
tion of 0.26 and 0.05 mag in the F336W and F110W bands, 
respectively. The former upper limit corresponds to an ab- 
solute brightness of M3oo,ab > —18.3 mag in the 300-nm 
rest frame, assuming F„ oc v^ 1 ' 2 at that wavelength. The 
NIR observation in F110W can constrain the rest frame op- 
tical/NIR brightness of the host galaxy, mostly attributed 
to late type stars and the old stellar population. The up- 
per limit corresponds to M > —18.2 mag in the 780 nm, 
assuming F v oc f -1 ' 2 . As is the case in the UV, this 
host is amo ng the faintest hosts in the lum i nosity dis- 
tributions of Savaglio. Glazebrook fe Le Borgnd (J2009J) and 
lLaskar. Berger fe Charvl (|201lh . 



9 This estimate assumes a 
locations inside the host 
licity requirement for the 



random distribution of the GRB 
jjalaxy. However, the low metal- 
GRB progenitor formation (e.g.. 



iHirschi. Mevnet fc Ma cdcr 2005) would increase the probability 
of GRBs being located in the outskirts of their hosts, given the 
typical metallicity gradients. 



http : //swift . gsf c . nasa . gov/docs/swif t/archive/grb_table/ 
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6 SUMMARY AND CONCLUSIONS 

We collected all observational data of GRB 070125 available 
in the literature, in order to investigate the environmental 
properties of the burst, and in particular to test the putative 
halo origin of this burst based on the few weak absorption 
lines i n the spectrum and the absence of an underlying bright 
host (|Cenko et al.ll2008t ). 

The NIR-to-X-ray SED at 1.9, 2.5 and 2.9 days after the 
burst revealed a synchrotron cooling break in the UV. Along 
the line of sight in the host galaxy, little reddening {E(B — 
V) = 0.016 ±0.007) is caused by a LMC- or SMC-type dust, 
while a low equivalent hydrogen column density cannot be 
constrained better than iV(H)= 15.5ii|;g x 10 20 cm -2 . The 
NIR-to-X-ray SED results are inconsistent with broadband 
SED analyses including radio observations, indicating that 
the assumptions in modelling efforts for GRB 070125 so far 
have probably been too simplified. More complex afterglow 
models should be explored, in particular to unify the light 
curves across the broadband spectrum with the NIR-to-X- 
ray SED. 

The analysis of the high S/N, but rather featureless 
FORS spectrum showed weak Mg n and C IV in absorp- 
tion, providing log A(Mgll) = 12.96io.il- Tne constraint 
on the Mg II column density, together with the upper limits 
for the other ions and the evidence for low metallicity from 
the SED, suggested that the GRB absorber is most likely a 
0.03 Zq < Z < 1.3 Zq Lyman limit system, also observed 
along the line of sight to five other GRBs. The comparison 
of the GRB 070125 measureme nts with the rest-fra me EWs 
of the GRB afterglow sample of lFvnbo et alj IJ2009T ) revealed 
that GRB 070125 is not unique, but, simply falls at the low 
end of the EW distribution. Furthermore, we selected two 
sub-samples of afterglows possibly similar to GRB 070125, 
the spectral-line poor and the low iV(H). Most spectral-line 
poor afterglows show a low Ay and low iV(H), suggesting 
that a short distance in the host galaxy is traversed by the 
line of sight. In any case, a quite featureless spectrum is not 
as uncommon for GRB afterglows as previously claimed. 

We demonstrated that the few and weak features in the 
spectrum of GRB 070125 can be the result of the particularly 
intense afterglow radiation ionizing a low log A(H i) ~ 19 
through a line of sight traversing not more than about 0.5 
kpc of its faint host galaxy. We showed that few and weak 
spectral features are not so uncommon in GRB afterglows, 
including the nearby GRB 101219B at z = 0.55, with asso- 
ciated SN, for which a sub-luminous host galaxy has been 
claimed. Finally, the non-detection of a host galaxy at the 
GRB 070125 position indicates that an underlying galaxy, if 
present, must be among the faintest GRB host galaxies. 

Given that about half of the GRBs in the spectral-line 
poor sub-sample have a detected host galaxy, it seems likely 
that these bursts are also located in the outskirts of a gas- 
rich, massive star-forming region inside its small and faint 
host galaxies, rather than in a halo environment. 
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APPENDIX A: 



GRB 070125 and the spectral-line poor absorbers 17 

Table Al. Swift UVOT photometry of the afterglow, corrected for Galactic extinction. 
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Table A2. NIR-to- X-ray SED modelling with a single power law (PL), a broken power law (BPL) and a t ied broken power law (TBPL), 
where the spectral slopes are tied to differ by A/3 = 0.5, as expected for a synchrotron cooling break GSari. Piran fc Narayanl ll998), 
Individual epochs (1.0, 1.4, 1.9, 2.5 and 2.9 days after the trigger) are modelled separately. Quoted errors refer to a 90% confidence level. 
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